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ABSTRACT 



We report the discovery with the Atacama Pathfinder Experiment (APEX) of an outflow of molecular gas in the radio-loud Seyfert 
galaxy IC 5063 (z = 0.01 10). In addition to the emission of the large-scale CO disk, a prominent blueshifted wing is observed in the 
C0(2-l) spectrum. IC 5063 represents one of the best cases of a fast jet-driven HI (and ionized gas) outflow, which is located at the 
site of a radio-bright feature about 0.5 kpc from the nucleus. It is possible that the blueshifted part of the molecular gas is associated 
with this outflow and is accelerated by the interaction with the radio jet. The outflow of molecular gas is characterized by an H2 mass 
of the outflowing component of between 2.25 ± 0.70 x 10^ Mq and 1 .29 ± 0.40 x 10* Mq and a mass outflow rate between 22 and 129 
Mq yr"' depending on the assumption for ax and assuming a luminosity ratio L'CO(2-1)/L'CO(1-0) = 1. This confirms that this may 
indeed be the dominant component in outflows driven by the nuclear activity that are also found in other objects. However, this high 
mass outflow rate cannot easily be supported for a long time, suggesting that the gas outflow in IC 5063 happens in bursts and is in a 
particularly strong phase at present. Owing to its proximity, IC 5063 serves as an excellent laboratory for understanding the impact of 
radio jets on the gas-rich inter-stellar medium. 
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1. Introduction 

The mechanism of energy feedback is considered to be the key 
to successfully model galaxy evolution (Silk & Rees 1998, di 
Matteo et al. 2005). In early-type galaxies, feedback is mainly 
provided by the energy released by the active nucleus, and gas 
outflows are one of the main signatures of this process in action 
(see Fabian 2012). 

Outflows driven by nuclear activity can be initially launched 
from the accretion disk or dusty torus that surrounds the black 
hole (BH), in the form of radiatively driven winds. The other 
main mechanism often considered is acceleration of the gas via 
radio jets/lobes. Understanding the relative importance of these 
mechanisms is the goal of many recent studies. Outflows have 
so far mainly been traced using ionized gas (see e.g. Crenshaw, 
Kraemer & George 2003, Nesvadba et al. 2007, Holt et al. 2009, 
Reeves et al. 2009, Tombesi et al. 2012, Harrison et al. 2012). 
Outflow effects have often been found to be confined to the in- 
nermost nuclear regions with relatively modest associated mass 
outflow rates (see e.g. Holt et al. 2011). Although this could 
have been a problem with the initial models of radiatively driven 
outflows (di Matteo et al. 2005), the possibility of a more com- 
plex, two-phase action on the inter-stellar medium (ISM) clouds 
makes this mechanism very relevant in the feedback process 
(Hopkins & Elvis 2010). 

Interestingly, more massive outflows have recently been dis- 
covered from the study of atomic neutral H I and molecular gas 
(e.g. Morganti et al. 2005a, 2005b, 2010; Femglio et al. 2010, 
Alatalo et al. 2011, Dasyra & Combes 2011, 2012). In a num- 
ber of objects, radio plasma jets have been suggested to play 



the dominant role in driving the outflows. Indeed, in some cases 
(Morganti et al. 2005b, Oosterioo et al. 2000), the location of 
the outflow appears to be clearly co-spatial with prominent radio 
features (e.g. hot spots). This has also been independently con- 
firmed by observations of the ionized gas component. Radio jets 
can indeed provide a particularly suitable and fast way of trans- 
porting the energy because they couple efficiently to the ISM and 
produce fast outflows from the central regions as required from 
feedback models (Wagner & Bicknell 2011, Wagner, Bicknell, 
Umemura 2012). It is surprising, however, that despite the en- 
ergetic interaction, these off-nuclear outflows can still have a 
component of relatively cold gas (< 1000 K) as detected in HI. 
Whether this is a common characteristic is an open question and 
is the motivation for the observations presented in this letter. If 
radio jets are as important - in certain conditions - as claimed by 
Wagner & Bicknell (201 1) and Wagner et al. (2012), the charac- 
teristics of the cold gas can provide constraints on the interaction 
physics. 

Here we present the results of CO(2-l) observations of 
the nearby, southern radio-loud Seyfert galaxy IC 5063 (z = 
0.01 lOfl selected because it represents one of the best examples 
of fast (~ 700 kms ') outflow of neutral hydrogen and ionized 
gas (Morganti et al. 1998, 2007, see Fig. 1 for an overview) that 
is located at about 0.5 kpc from the nucleus and is coincident 
with the bright radio lobe (Oosterioo et al. 2000). Therefore, the 
mechanical energy associated with the radio jet/lobe has been 
suggested to be the main mechanism that drives the outflow. This 
object is, therefore, an ideal candidate for understanding the de- 
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' Throughout this paper we use the Hubble constant Ha= 70 km s"' 
Mpc"' and = 0.7 and £1m = 0.3. At a distance of IC 5063, this 
results in 1 arcsec = 200 pc, (Wright 2006). 
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Fig. 1. Some relevant characteristics of IC 5063 (taken from 
Morganti et al. 1998). Left: the morphology of the radio contin- 
uum emission (white contours) is superimposed on an HST im- 
age. The contours show the bright radio lobe (the W component), 
the core, and the fainter E lobe. The APEX CO(2-l) observation 
is centered on the radio core but with a beam of 30", it covers the 
entire radio source. The right figure illustrates the H I position- 
velocity plot taken along the major axis of the H I disk (black 
dashed line on the radio image). The large-scale H I disk in emis- 
sion is shown with solid contours while the broad blueshifted 
component detected in absorption (dashed contours) is clearly 
evident, covering velocities well outside the range of the normal 
gas disk (Morganti et al. 1997). The spatial axis in this plot cov- 
ers the galaxy out to 2.5 arcmin radius. VLBI follow-up H I ob- 
servations have shown that the blueshifted absorption is located 
at the position of the western brighter radio lobe, about 0.5 kpc 
from the nucleus, suggesting that the jet-ISM interaction occurs 
at this location (Oosterloo et al. 2000). 

tails of the gas associated with such interactions and its relevance 
for feedback. 

IC 5063 was previously observed and detected with SEST in 
CO(l-O) by WikHnd, Combes & Henkel (1995). The availabihty 
of these data prompted us to apply for observations of the CO(2- 
1) component. The CO(l-O) profile reported by Wiklind et al. 
(1995) is, however, relatively narrow (AvfWHM - 163 kms"') 
and offset in velocity compared, e.g. to the H I emission profile. 
In early-type galaxies, molecular and atomic neutral gas tends 
to show similar kinematics (Davis et al. 2011). Thus, the new 
observations can also serve as a comparison for this puzzling 
difference. 

2. C0(2-1) in IC 5063: the APEX observations 

Observations with the Atacama Pathfinder Experiment (APEX) 
12-m antenna were conducted in two campaigns in 2008 us- 
ing the APEX-1 instrument with the FFTSl backend tuned to 
230.538 GHz, the frequency of the CO(2-l) Une at the red- 
shift of IC 5063 (z = 0.01 10). The central position of IC 5063 
(RA=20h52m02.3s Dec=-57d04m08s) was observed on August 
8-9, 2008 in good weather conditions (for these frequencies). 
The precipitable water vapor (PWV) was < 2 mm. IC 5063 was 
observed for a total integration time of 146 minutes. 

The observations were made using 4096 channels, resulting 
in a velocity resolution of 0.32 km/s. The two fast Fourier trans- 
form Spectrometer (EFTS) units almost fully overlap, resulting 



in a total bandwidth of 1 GHz (~ 1250 kms '). For the individ- 
ual scans of each of the EFTS units, a linear baseline was sub- 
tracted before adding all spectra. The final spectrum is smoothed 
to ~ 15 kms"' bins. The rms-noise at 15 kms"' resolution is 
~ 1 mK. The data were reduced with the CLASS software from 
the Gildas package (|http://www.iram.fr/IRAMFR/GILDAS|). 

At the frequency of these observations, the spatial resolution 
of APEX is ~ 30" (i.e., covering a radius to ~ 3 kpc from the 
center). At this resolution, the central beam includes the whole 
radio continuum source and the inner part of the gas disk ob- 
served in H I emission (see Fig. 1 for reference). Inside this ra- 
dius, the H I emission covers the full velocity range of the rotat- 
ing H I disk (because of the steep rotation curve typical of early- 
type galaxies) and also includes the broad, blueshifted compo- 
nent observed in HI absorption against the radio source (see 
FiglD. 

The final CO(2-l) profile observed in IC 5063 is shown in 
Fig.|2] The new result that is immediately apparent from this fig- 
ure is the large width of the CO(2-l) profile: about 800 km/s 
FWZI. Thus, the CO(2-l) profile appears to be significantly 
broader than what was derived in CO(l-O) by Wiklind et al. 
(1995). The width and the central velocity of the CO(2-l) Hne 
(3400 vs 3349 km/s in CO(l-O)) appears to be more consistent 
with the HI profile (as can be seen in Fig. 2). Furthermore, the 
CO(l-O) does not appear to be asymmetric and does not have a 
blueshifted wing. 

It is perhaps more interesting that the profile is also much 
broader than the integrated H I emission profile obtained from 
the entire H I disk. This is illustrated in Fig. 2 where the H I emis- 
sion profile is taken from Morganti et al. (1998) and represents 
the full velocity amplitude of the HI gas in the large-scale disk, 
which is unaffected by the active galactic nucleus (AGN). Also 
indicated in Fig. 2 is the velocity range corresponding to the 
H I outflow seen in absorption. The important result is that the 
CO(2-l) spectrum is asymmetric with a clear blueshifted wing 
that extends well outside the velocity range of the gas disk and 
into the velocity range of the H I outflow. The fact that the CO(2- 
1) profile is broader than the H I emission profile illustrates that 
an extra component must be present, with velocities exceeding 
those of the gravitationally rotating gas. This suggests that the 
molecular gas is participating in the same outflow that drives out 
the HI. 

The bandwidth of the CO(2-l) observations to some extent 
limits the sampling of the full range in velocities covered by the 
H I absorption. Therefore, the full width of the blueshifted com- 
ponent of the molecular gas is somewhat uncertain. However, the 
data clearly show that excess emission that cannot be explained 
by poor subtraction of the linear baseline during the reduction 
steps. 

3. A new component of the gaseous outflow: the 
molecular gas 

The blueshifted wing in the CO(2-l) profile suggests that in 
IC 5063 we are observing the molecular gas associated with 
an outflow in addition to the molecular gas associated with the 
large gas disk. Despite the limitations present in the CO(2-l) 
data mentioned above, we now derive some parameters to char- 
acterize the outflow of molecular gas. 

We estimated the flux of the outflow component by fitting 
two Gaussians to the CO profile. The result of the fitting is shown 
in Fig. 3. In this way we obtained Ico - 0.54 ± 0.17 K kms"'. 
We used a conversion factor of 39 Jy K"' for APEX. To es- 
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Fig. 2. Profile obtained from the APEX CO (2-1) observations. 
The profile is extremely broad (almost 800 km s"' FWZI) and 
asymmetric, showing a pronounced blueshifted wing. The pro- 
file reaches velocities well outside the range covered by the regu- 
larly rotating large-scale gas disk. This is illustrated by the com- 
parison with the velocities covered by the H I emission (indicated 
in gray in the plot, with arbitrary intensity scale). The blueward 
wing of the CO(2-l)) profile covers velocities similar to those of 
the H I outflow detected in absorption. 



of ~ 22 Mq yr" (for the conservative assumption of ax - 0.8 
Mq), reaching up to 129 Mq yr"' using the result of the Gaussian 
fitting and ax - 4.6 Mq. 

For these relatively high mass outflow rates, one may won- 
der whether they can be easily supported in the case of IC 5063. 
The total mass of the molecular gas - in the central beam - that 
follows the regular rotation is 4.96 + 0.50 x 10** Mq (based on the 
results of the Gaussian fitting and using ax - 4.6 M©). Of this, 
we can assume that about a tenth (~ 5 x 10^ Mq) is included in 
the central kpc. This is derived by very roughly assuming that it 
follows the ratio between the beam area (radius 3 kpc) and the 
area affected by the jet (radius ~ 1 kpc). Due to the collimated 
nature of the outflow, only part of the gas in the inner region is 
affected by the jet at any instant. However, because of the rota- 
tion of the inner gas disk, eventually the entire gas of the inner 
disk will be affected. The orbital timescale in the inner region is 
about 3x10^ years, implying an average mass outflow rate of 
~ 1 Mq yr"'. This is, of course, a very simplified view, and the 
mass of the molecular gas in the central beam could be higher 
if an exponential disk distribution were assumed. However, even 
for a substantially larger portion of molecular gas, the mass out- 
flow rate would stay below what was measured. This difference 
may suggest that the gas outflow in IC 5063 occurs in bursts and, 
at present, is in a particularly strong phase. If this is the case, this 
would have interesting implications for an explanation of some 
of the peculiar characteristics of this source, and we return to 
this question in Sec. 5. 



timate the H2 masses we assumed a luminosity ratio L'CO(2- 
1)/L'CO(1-0) - 1 typical of thermalized and optically thick gas, 
as could be the situation in denser and hotter gas associated with 
an outflow. If this is the case, it could explain why the broad 
component was not detected in Wiklind et al. (1995). However, 
the physical conditions of these outflows are still poorly known 
and other possibilities have been proposed (Alatalo et al. 2011). 

The choice of the conversion factor ax - Mh^ / L'^^ Mq (K 
kms 'pc^)"' is not easy because very little is known about out- 
flows (see discussion in Alatalo et al. 201 1). To be conservative, 
we used the two extreme values of ax - 0.8 Mq (K km s^'pc-)"' 
derived by Downes & Solomon (1998) for ultra-luminous in- 
frared galaxies and the standard ax - 4.6 Mq (K kms 'pc^)"' 
(see Solomon & Barrett 1991 for a discussion). The derived mass 
of the molecular outflow ranges from 2.25 + 0.70 x 10^ Mq to 
1.29 + 0.40 X 10'^ Mq. Thus, the derived mass of the molecu- 
lar outflow appears to be higher than that deduced for the H I 
outflow (i.e., 3.6 x lO*" Mq; Morganti et al. 2007). 

The total mass of the molecular gas - in the central beam 
- following the regular rotation (i.e., centered on the systemic 
velocity) and estimated using similar assumptions as described 
above, is 4.96 + 0.50 x 10^ Mq. This means that the ratio be- 
tween the molecular masses derived from the CO(2-l) and the 
one quoted by Wiklind et al. (1995) and derived from the CO(l- 
0) is 1.3. 

Following Alatalo et al. (201 1), we derived the mass outflow 
rate of the molecular gas. The spatial resolution of the APEX 
observations do not allow one to identify the exact location of 
the molecular outflow in IC 5063. If the molecular outflow is 
located in/near the western and radio-bright lobe, i.e., about 0.5 
kpc away from the nucleus, the time for the outflow to reach this 
radius is t^„, ~ 1 Myr, assuming the highest velocity observed 
for the CO (~ 400 km s"' relative to the systemic velocity). The 
mass outflow rate estimated as M = M/t^„, has a lower limit 



4. Origin of the outflow 

Existing data on the observed outflow of ionized and atomic gas 
in IC 5063 suggest that these outflows are driven by the interac- 
tion between the radio jet and the ISM. The molecular compo- 
nent could be part of the same outflow. 

Molecular gas associated with the jet/ISM interaction has 
also been suggested on the basis of observations of warm molec- 
ular gas by NICMOS (Kulkarni et al. 1998) and, more recently, 
in the mid-IR by Spitzer (Guillard et al. 2012). KuUcami et al. 
(1998) found that the [Fell] /II. 644 fim and H2 i2.1218 yum 
emission presents a one-to-one spatial correspondence with the 
radio lobes, suggesting that these lines originate in shocks pro- 
duced by the advancing radio jets. Furthermore, the asymme- 
try in the [FeII]/H2 ratio between the eastern and western sides 
suggests an excess of molecular gas on the W side, where the 
radio lobe is brighter This reinforces the idea of the jet interact- 
ing with a particularly rich and dense medium. More recently, 
this scenario has obtained support from Spitzer mid-IR data. 
The warm H2 appears to be particularly bright in radio galax- 
ies with an HI outflow, and most likely the dissipation of the 
kinetic energy of the radio jet is powering this bright H2 emis- 
sion. IC 5063 is one of such cases (Guillard et al. 2012). As de- 
scribed in Guillard et al. (2012), and following theoretical stud- 
ies (Meflema et al. 2002, Gaibler et al. 2012, Wagner & Bicknell 
2011, Wagner et al. 2012), the jet itself cannot affect a signifi- 
cant fraction of the ISM of the host galaxy because it is very col- 
limated. Instead, the expanding cocoon of hot and tenuous gas 
inflated by the jet (e.g. Begelman & Cioffi 1989) could trans- 
fer part of its kinetic energy to the surrounding ISM by driv- 
ing shocks and/or by turbulent mixing. Numerical simulations 
have studied the effects of a cocoon of shocked material run- 
ning over a clumpy medium and/or a massive molecular cloud 
(Mellema et al. 2002, Gaibler et al. 2012, Wagner & Bicknell 
2011, Wagner et al. 2012). A compression phase of this pre- 
existing gas is expected, followed by fragmentation and cooling. 
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Fig. 3. Fit of the CO(2-l) profile showing two Gaussian compo- 
nents. See text for details. 

This process generates dense, cool, fragmented structures that 
can be entrained and accelerated by the cocoon that surrounds 
the radio jets/lobes. The key to these models is that radiative 
cooling must be very efficient to be able to produce accelerated 
clouds of cold gas (which would be observable in H I and CO). 
In this scenario, H I would represent an intermediate phase in the 
cooling process of the warm gas, and the molecular gas would 
be the final stage. The simulations also showed that the transfer 
efficiency of kinetic energy and momentum from the jet to the 
dense {hh > 10^ cm"-') gas can be high (10-70%). 

In summary, the detection of cold gas and the expectations 
of the simulation underline the complex and multi-phase nature 
of the outflows. 

5. Discussion and conclusions 

We have presented the detection of molecular gas outflow in the 
radio-loud Seyfert galaxy IC 5063. This represents another case - 
in addition to Mrk 23 1 (Feruglio et al. 201 0), NGC 1 266 ( Alatalo 
et al. 201 1), and 4C12.50 (Dasyra & Combes 2012) - where fast 
outflows of molecular, atomic neutral and ionized gas are found 
to co-exist in the same object. 

In addition to the possibility that these outflows are driven by 
ffie interaction between the radio jet plasma and the rich ISM, it 
is also interesting to consider the possibility that the same in- 
teraction has a major impact also on the characteristics of the 
radio source and, in particular, on enhancing the radio emission. 
This could also explain some of the intriguing characteristics of 
IC 5063 described above (i.e., high radio power for a Seyfert 
galaxy, large asymmetry between the two radio lobes, high mass 
outflow rate at the present stage). Because of the compression 
of the magnetic field and the increased density of particles, the 
radio luminosity can be boosted by the interaction (Tadhunter 
et al. 2011). This has been suggested for the young, compact 
source PKS 1814-637 hosted by a disk galaxy with a rich ISM 
(i.e., H2, and PAH emission along with HI and silicate absorp- 
tion features) — characteristics that are quite rare among pow- 
erful radio sources. Because of this, PKS 1814-637 and other 
similar objects have been suggested by Morganti et al. (201 1) to 
be something like impostors: objects intrinsically of low power 
that are selected in radio-flux-limited samples because of the ef- 
ficient conversion of jet power into radio emission. These are 
objects with a rich, dense ISM that are expected to have unusu- 
ally strong interactions between the jets and the ISM (Tadhunter 
et al. 201 1). They may also represent a missing link between ra- 
dio galaxies and radio-loud Seyfert galaxies. IC 5063 is one of 
the most radio-loud Seyfert galaxies. Thus, it is intriguing that 
in this object we see clear evidence that the radio jet strongly in- 



teracts with the ISM, in agreement with what is expected in the 
above scenario. The suggestion that the mass outflow rate is, at 
present, particularly high in IC 5063, additionally supports that 
ffiis object is in a special phase at the moment. 

Finally, although it is massive, the cold gas component of the 
AGN-driven outflows does not always appear to be large enough 
to match what is required by quasar-feedback models (di Matteo 
et al. 2005). In IC 5063, the kinetic power derived by combining 
the H I and CO components (assuming they are tracing the same 
outflow) is between 7.8 x 10'*^ and 1 .7 x 10'*-' erg s"\ depending 
on the mass outflow rate assumed for the molecular gas. These 
estimates correspond to ~ 10 - 20 % of the nuclear AGN bolo- 
metric luminosity (as derived in Morganti et al. 2007), which 
suggests that the impact of the outflow on the ISM is relevant in 
the evolution of the host galaxy. However, one should consider 
the possibility that the mass outflow rate is particularly high at 
present but cannot be maintained at this level for a long time, 
and therefore may have an effect only for a limited period in the 
life of the galaxy. 

The next step in the study of this object is to use higher reso- 
lution and more sensitive observations to find the exact location 
with respect to the other components of the molecular outflow, 
and investigate its characteristics in greater detail. In this way, it 
will be possible to verify the scenario sketched here. 
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